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Abstract 

We discuss one-particle inclusive B decays in the limit of heavy b and c 
quarks. Using the large-A'^c limit we factorize the nonleptonic matrix 
elements, and we employ a short distance expansion. Modeling the 
remaining nonperturbative matrix elements we obtain predictions for 
various decay channels and compare them with existing data. 

PACS numbers: 13.25.Hw, 12.39.Hg 



1 Introduction 



Over the last ten years, methods have been developed to describe the decays 
of heavy hadrons within a framework based on the heavy mass expansion of 
QCD. Using this approach, model dependences have been drastically reduced. 

The heavy mass expansion has been formulated for exclusive as well as 
for inclusive decays. While in the former case the heavy mass expansion 
is the well-known formalism of heavy quark effective theory (HQET) 
the latter case additionally requires a short distance expansion (SDE) 
which is very similar to the operator product expansion (OPE) used for deep 
inelastic scattering. Both approaches are QCD based methods and have a 
good theoretical foundation. 

However, nothing comparable exists for for a theoretical description of 
one-particle inclusive decays such as i? — > DX or B ^ R'X. Unlike the fully 
inclusive case, the expression for the rate always involves a projection on a 
specific particle, thereby spoiling a straightforward short distance expansion. 
The semileptonic case has been discussed in 0, and the present paper is 
devoted to a study of the nonleptonic decays. 

We shall concentrate on decays of the form B DX and B DX. 
These decays have already been considered long ago in the context of models 
01 . Here we aim at a QCD description and exploit the heavy mass limit for 
both the b and the c quark. Furthermore we shall use the large- iV^ limit to 
factorize the hadronic matrix element. In this way we can identify the parts 
which allow a short distance expansion. The remaining contributions have to 
be parametrized and we shall discuss simple forms of this parametrization. 

We first discuss the right charm contribution arising from quark level 
b — » cX decays and the wrong charm contribution arising from b cX 
separately. After recalling some facts about the phenomenology of D mesons 
in Sec. ^, we predict rates and spectra for various decay channels in Sec. |^. 
Finally we compare the results with existing data. 

2 The effective Hamiltonian 

and the right charm contribution 

The relevant effective Hamiltonian for the decays B DX is given by 

H.ff = Hi';] + Htll (1) 
where the semileptonic and nonleptonic pieces are 

H'fl] = % {bc)v-A Mv-A + H.C., (2) 
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y^q=d,s\ k=l k=l k=3 / 

(3) 

Oi,2 and Pi,2 are the current- current operators, and O3...6 are the QCD pen- 
guin operators. For a full list of operators, see, e.g., the review We shall 
only consider the Cabibbo-favored decays and neglect penguin contributions. 
Hence the operators we need are 



= Oi = ibc)v-Aiud)v-A 

(4) 



= 02 = {bT''c)v-Aiun)v- 



= Pi = {bc)v^A{cs)v^A 
= P, = {hT'^c)y^A{cT'^s)v-A. 

Ci and C2 are the Wilson coefficients encoding the short distance physics 
and are the generators of color-SU(3). The operators Oi and O2 as well 
as the semileptonic Hamiltonian H^J'f] contribute to right charm transitions 
only, while Pi and P2 contribute to both right and wrong charm processes. 

It is well known that in the large-A^c" limit the matrix elements of the 
four-fermion operators factorize into products of two current matrix elements 
0. Factorization has also been investigated on a phenomenological basis 
and found to work well for exclusive nonleptonic decays of B mesons [|^. 
The contributions of O2 and P2 vanish in the factorization limit since the 
currents are color octets. Furthermore, the coefficient Ci differs from unity 
only through radiative corrections which are very small and will be neglected. 

The semileptonic case has already been studied in so we focus on the 
nonleptonic modes. We consider the matrix element 



G(M2) = E {D^'\pj,)X\H^jf\B{pB)) 

X 
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{2^Y6\pb-Pd-Px) (5) 



where the states \X) form a complete set of momentum eigenstates with 
momentum px and -f^e// is the relevant part of the weak Hamiltonian. The 
function G depends on the invariant mass 

M' = {ps-pj^f (6) 

of the state \X). It is related to the decay rate under consideration by 

dT{B D^*^X) = d^jy G{M^), (7) 

U\ 

where d^jy is the phase space element of the final state D meson. 
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Due to the different final states cud and ccs, there are no interference 
terms between Oi and Pi. The contribution of the channel b cud is 

Gi(M^) = ^ \V^,V:f j:i^7rmM-Px) (8) 
^ X 

{B{pn)m,{l - ^,)b)(dY{l - i,)u)0*\pT^)X) 
(D«(p^)X|(n7^(l -75)rf)(H(l -75)c)|i?(pi.)). 

Using the factorization of the large- iVc limit we have (see Fig. |1|) 

Gi{M^) = ^\Vciy:f\Ci\' EE(2vr)V(M-px-px') (9) 

^ XX' 

{B{pB)\{ci,{l-l,)h)\D''*\pD)X) mW-l^)u)\X') 
{X'\{uY{l - 75)rf)|0) {D''*\pj^)X\{hu{l - 1,)c)\B{pb)). 

It is convenient to define two tensors 

K,,{PB. M, Q) = Y.{2'kY5\M -Q-Px) (10) 

X 

(S(Pb)|(c7,(1 - ^,)b)\D^*\p^)X){D^*\p^)X\{hua - 75)c)|5(pb)) 

and 

P,AQ) = E(2vr)V(Q-px')(0|(37;.(l-75)^)|X')(^'l(^7.(l-75)^^)|0) (11) 

X' 

in terms of which we obtain for the rate 

Gi{M') = ^\Vc,V:f\Ci\'J -^^K,Apb,M,Q)P^^{Q). (12) 
The tensor P^j, involves only light quarks and can be rewritten as 

P,.iQ) = I rf'a;e-^«^(0|(d(x)7^(l - i,Mx))iuiO)l.i^ ' 75)ci(0))|0). (13) 

For sufficiently large Q this quantity has a short distance expansion in inverse 
powers of Q. However, the momentum Q is not a measurable kinematical 
quantity. In particular, it is not equal to the recoil M , but in most of the 
phase space, it should be of the same order as M, namely 0{mb — rric)- Close 
to the nonrecoil point, M is large, corresponding to either large Q or a large 
momentum of the gluon depicted in Fig. |I[ While in the former case the OPE 
treatment is justified, the latter case could be treated perturbatively. In other 
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Figure 1: Momenta and velocities used for the description of the right charm 
contribution. 

words, the two expansion parameters for the total rate are I^qcd/ (^b — '^c) 
and as{mc)/TT. However we also consider the momentum spectra of the D 
mesons. In that case the expansion parameter depends on the energy of the 
jD meson since 

Q'^ pa M'^ ^{m^-mcf -2m^{Ec-mc). (14) 
The leading term of the short distance expansion yields the partonic result 

PAQ) = ^ {Qi^Qr^ - 9,^^Q^) e»(Q') (15) 

where we have assumed both light quarks to be massless. 

For the quantity K^i, it is convenient to use the heavy mass limit for both 
the bottom and the charm quark. To this end we redefine the phases of the 
quark fields as 

h{x) = K{x) e-^"'*'^^ c{x) = c,/(x) e-^"^=^'^ (16) 

where the velocities are defined as pb = tubv and = m-^' . In the 
following we shall work in the infinite mass limit for the h and the c quark, 
and we obtain 

K^,{m,v,M,Q) = I ^expH(M - Q)z\ (17) 

{B{v)\{cAz)l,{l-l,)K{z))\D^*\v')X) 
0*\v')X\iMOhA^ - 7,)c.'m\B{v)), 
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where in the heavy mass hmit M becomes — irtcv' . 

Using the result in Eq. ( pTSf ) for P^j^ we can express the rate in terms of 
the two quantities 



K,{v,v',Q)=K^^^{mbV,M,Q), 




(18) 



(19) 



and obtain 




(20) 




jQ' Qm[K,{v,v',Q) - K2{v,v',Q)] 



Not much can be said about the functions Ki and K2. They do not 
have an obvious short distance expansion due to the projection on the final 
state D meson. The only restriction we have is from the spin symmetry of 
the heavy b and c quarks. To implement these symmetries we use the spin 
projection matrices for the heavy mesons 



In the heavy mass limit, D and D* become degenerate and constitute the 
ground state spin symmetry doublet of the D meson system. For phenomeno- 
logical applications we shall later take the mass splitting between the D and 
D* mesons into account, although formally this is a l/rric effect. 

As far as the spinor indices are concerned, the tensor K^^i^ is given by 



where the remaining indices are light quark indices which have to be con- 
tracted using the most general four- index object. 

Since there are many possibilities to contract the indices, the discussion of 
the general case would leave us with a large number of unknown functions, 
so we have to make a choice. The matrix elements appearing in K^^iy are 
identical to the semileptonic case. Therefore we shall use the same ansatz 
for the nonleptonic as for the semileptonic case 0] and write 



Hb{v) = -7m^(l + 1^)75, 



(21) 



%*)K) = ^v^(l+^') 



75 (pseudoscalar meson), 
/ (vector meson). 



K^^ipB, M, Q) = {2tiY5\M - Q) ri{vv') 
Tr[:ffB(i;)7^(l - 75)i^75w(i^')] Tr[:ff^M(t;07.(l - 75)i^B(t^)]. 



(23) 
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In this way we can model the two functions Ki and K2 in terms of the single 

nonperturbative function r]{vv'), the advantage being that spin symmetry 
relates the rates of i? ^ DX and B ^ D X to each other. 

Finally we also have to consider the right charm contribution of the op- 
erator Pi, which is 

G^iM') = ^ \V^V^\' \C^\' Y.i^nrS\M-px) (24) 
^ X 

{B{ps)m,ii - 7^msrii - i,)c)\d^*\p^)x) 

{D^*\p^)X\{cY{l - l,)s){hu{l - 1,)c)\B{pb)) . 

The calculation is exactly the same as before, but the short distance ex- 
pansion of P^i, yields a different result since now a massive charm quark is 
involved. We have 

P',AQ) = E(2vr)^5^(g-pxO(0|(57,(l-75)c)|X')(X'|(c7.(l-75)^)|0) (25) 

X' 

from which we have the leading order contribution 

p;,(g) = (A(g2) q^q^ _ b{q^) g v) ©(g' - (26) 

with 



Bm = 




and we get 

G2{M^) = ^ \V.,V:f \C,\' J e(g^ - ml) (28) 

[A{Q^)K,{tKv\Q) B{Q')K2{v,v',Q)]. 
Thus the same two nonperturbative functions appear in G2- 

3 The wrong charm contribution 

Wrong charm decays can only be mediated by the operators Pi and P2. In 
the large- A^c limit, the contribution of P2 can be neglected and we have 

Gs{pb,Pd) = ^ \V^V:f |Cip J2i^nr^'iPB - PD - Px) (29) 
^ X 

{B{pb)\{ci,(1 - 75)&)(s7^(l - 1,)c)\D^*\pd)X) 

{D^*\pn)X\{cY{l-l,)s){hA^-l^)c)\B{pB)). 

6 




Figure 2: Momenta and velocities used for the description of the wrong charm 
contribution. 



Again using the factorization of the large- A^^c limit we get (see Fig. |^) 
G;{pb,Pd) = ^ \V^V:f \C,\' Y^Y.C^rrr^'iPB-PD-Px ~Px') (30) 

^ XX' 

{B{pB)m,{i-i,)h)\x) {x\{Tnu{i-i,)c)\B{pB)) 

(0|(^7'^(1 - 75)c)|I)W(Pd)X') {D^*\pn)X'\{cY{l - 75)s)|0). 
In a similar way as before it is convenient to define two tensors 

Q) = E(2^)'^'(PB -Px-Q) (31) 

X 

{BiPB)m,a-i5mx){x\{hu{i-i5)c)\BiPB)) 

and 

R,u{pD,Q) =Y.i'^7^y^'iQ-PD-PX') (32) 
X' 

(01(^7.(1 -75)c)|/^W(Pd)X')(Z}W(pz,)X'|(c7.(1-75)s)|0), 
in terms of which the rate becomes 

GsipB,PD) = ^\VctV:f\C^\' J ^,K'^Apb.Q)R'''{pd.Q). (33) 

The quantity K'^^ is fully inclusive and one may perform a short distance 
expansion. In the heavy mass limit for the h quark it is convenient to rescale 
the h quark field as in Eq. ([I6|) , and we obtain 

K'^Apb, Q) = J exp[-i(mbt; - Q)y] (34) 
{B{v)\{c{y)^,{l - 75)6.(2/))(&.(0)7.(l - 75)c(0))|i?(t;)). 
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In the region where the momentum rribV — Q is large, we can perform a short 
distance expansion. As before, m^v — Q is not an observable, but it is of 
the order of mf, — in most of the phase space. The leading term is the 
dimension three operator b^bi,, the matrix elements of which are normalized 
due to heavy quark symmetry. Thus we obtain 

K^MbV, Q) = 2Trd ({rribv - Qf - m^) (35) 

The other factor R^,, involves a projection on a D^*) meson in the interme- 
diate state and has to be parametrized. Heavy quark spin symmetry for the 
c quark implies that R^i, is of the form 

Rn„{mcv', Q) oc Hjj(»){v') 7^(1 - 75) O 7^(1 - 75) if^w(t;'), (36) 

where the remaining light quark indices have to be contracted using the most 
general four-index object. 

Again there is quite a large number of possibilities to contract the indices, 
making it useless to discuss the general case. Hence we shall only give two 
physically motivated ways for modeling this quantity. 

The first model ansatz corresponds to factorization: 

R^Am,v',Q)^2F{v'Q,Q^) (37) 
Tr{Hj,,.,{v') 7^(1 - 75)] TV[7,(1 - 75) Hj,,.,{v% 

the second one is inspired by the parton model and is defined through 

RA^cv\ Q)=2 F{v'Q, Q2) (38) 
Tr[^^„(^;')7^(l - 75)(^ - m,fh,{l - 'y,)Hj,M{v')]. 

We shall discuss the functions F and F using data. 

4 Effects of D* ^ D decays 

The right charm decays of 6 = +1 mesons are the ones into or D~ 

B+ ^ D^*)'^ X, 
B+ £)(*)- X, S° 

while the wrong charm decays are 

B+ X, 
B+ £)(*)+ X, S° 



D^*^' X, 
L>W- X, 



n(.*)o X, 

£)(*)+ X. 



(39) 



(40) 
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For the semileptonic case the charge of the the lepton tags the b flavor of the 
decaying B meson. Wrong charm semileptonic decays are suppressed by the 
large charm mass and will be ignored in our discussion. 

Off the heavy mass limit the degeneracy between D and D* mesons is 
removed. In fact the mass difference is large enough to allow strong decays 
into pions. Thus the rate for decays into D mesons is the sum of a direct 
contribution and the contribution arising from the decay chain B —>■ D*X' 
DX. In the narrow width approximation, the latter is obtained by weighting 
the rate for B —>■ D*X by the branching ratios D* DY where Y is either 
a pion or a photon. 

While D*^ decays are governed by the isospin Clebsch-Gordan coefficients 
receiving only tiny corrections from phase space effects and from the radiative 
process, 

Br(D*+ ^ D+Y) = 0.32, Br(D*+ ^ D^Y) = 0.68 (41) 

in Z)*° decays isospin invariance is maximally broken by phase space effects 
m such that 

Br(L)*° ^ D^Y) = 1. (42) 
Consequently the total branching ratios to pseudoscalar D mesons are 

Br (5 ^D-X) = Bi{B D^^^X) + 0.32 Bi{B D*-X) (43) 
Br(5^;D°X) = Br(5^D°i^X) + Br(fi^D*°X) +0.68Br(fi^D*-X) 

where D^ir refers to the direct contribution, i.e., the contribution where no 
D meson appears in an intermediate state. 

We shall later consider spin and charge counting in one-particle inclusive 
B decays and define the following ratios: 

Br{B-^'L)*X) 

rs = — ^ 44 

Br(5 ^ Ddir^) 

Bt{B^&X) 

^« = BriB^D-Xy ^''^ 

We assume that D* decay is the only relevant isospin violating effect. Under 
this assumption the charge counting ratio is governed by the spin counting 
ratio: 

l + 0.32r5' ^ ' 

Assuming equal rates for the pseudoscalar mesons and for each polarization 
state of the vector mesons, the naive expectations for these two ratios are 
rs = 3 and rg = 3. 
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Apart from the ground states D and D*, exited D mesons may be pro- 
duced in B meson decays. Even in the heavy mass hmit a finite sphtting to 
the ground state mesons remains, which means that these states will decay 
into D or D*. In order to analyze decays into higher resonances, one would 
need to study the state \X), which we take fully inclusive. Thus we may rely 
on parton hadron duality to obtain a correct description of the D and D* 
production up to leading order in the heavy mass expansion. 

Charmed mesons can also be produced via charmonium intermediate 
states. These states show up in the invariant mass distribution dV /drnj^j^ as 
more or less pronounced resonances. The distribution dT/dmcc is certainly 
different, but by global parton hadron duality and the assumption that all 
charm quarks hadronize into D mesons, their integrals are equal, so the char- 
monium contributions are included in our calculation. 



5 Results 

Even after restricting the number of possible form factors by the ansatze 
in Eqs. (p3D, (p^) and (|38D , we are still left with unknown nonperturbative 
functions. As far as the right charm contributions are concerned we shall 
first consider the decays into D mesons. According to Eq. (p3D, these are 
described by a single function rj{vv'). Spin symmetry relates these decays to 
the ones into D mesons, however, this is only true in the heavy quark limit, 
where again 7]{vv') is the nonperturbative input. In reality one has to take 
into account that D mesons decay into D's, and thus one has to add this 
contribution using Eqs. (^). 

Still all right charm decays are given in terms of the single function ri{vv'). 
Comparing the present case to the semileptonic one, we shall use the same 
saturation assumption as in P] and write 

^{vv') = \^{vv')\\ (47) 

For numerical calculations we use the measurement |^ of the Isgur-Wise 
function 

^{vv') = 1 - a{vv' - 1), a = 0.84. (48) 

Since the exclusive semileptonic decays are spectatorlike and since in the 
limit of factorization, the matrix elements of the heavy quark current are the 
same in the semi- and nonleptonic cases, dominance of spectatorlike decays is 
a natural assumption for the nonleptonic right charm well. Therefore 

we use the results of Sec. ^ for the channels 

^dir ^dir X, (AQ\ 

B+ X, D*- X, ^ ' 
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neglecting possible contributions from nonspectator channels such as 
D~X. Another class of decays allow to obtain a factorizable expression 
for ifeff after a Fierz transformation. However, the color indices need to 
be rearranged as well, yielding a suppression by one power of 1/Nc in the 
amplitude. Since we are working to leading order in 1/Nc, we find a vanishing 
rate for —>■ D*^X, whereas the channel B^ D^X is fed by the decay 
chain via the D*~ meson. 

For the contribution of the quark level decay b — >■ ccs, we assume the 
same number of right charm quarks to hadronize as D as wrong charm 
quarks hadronize as either D^*'' or D^*\ Using the CLEO wrong charm 
measurement Bt(B DX) = (7.9 ±2.2)% and the Particle Data Group 
average By{B DfX) = (10.0 ± 2.5)% [|lT| and neglecting a possible right 
charm B contribution, we obtain a total right charm contribution 

from b — >• CCS of about 18%. 

It is known that the channel b —>■ ccs receives large radiative corrections 



computed in [|I2[. We shall not include these corrections in their detailed 
form, rather we shall take into account their bulk effect by adjusting the 
charm quark mass in Eqs. (p6D-(p8|). Inserting an "effective" mass mf^ = 
1.0 GeV into the tree level relation, the measured rate is reproduced. 

Data are sparse for the wrong charm part, therefore we replace the un- 



known functions in Eqs. (^) and (|38D by 



F{v'Q,Q') = f2n5(iQ- 
F{v'Q,Q^) = f27i5hQ 



(50) 



where / and / are constants. This parametrization using a delta function is 
motivated by the negligible mass of the strange quark against which the two 
charmed quarks recoil. It can actually be checked by measuring the invariant 
mass uix' of the additional decay products in 5 — > DDX' which should turn 
out to be small. Setting / to one, Eq. ( pS]) reproduces the well-known parton 
model result. The fits yield 

/ = 0.147 GeV and / = 0.121. (51) 

Note that one naively expects / = 1/8 when assuming isospin invariance in 
B D*X and B Djir-^- Under this assumption, we have two charge 
states and four spin states, hence in total eight states contributing. 

For the semileptonic contributions, we use the results of ,[] but apply 
the same approximations as for the nonleptonic channels. We neglect the 

"'^Note the sign error afflicting one term in Eq. (43), which should read 
Ef^\y) = -%^MD*- -> ^X)\{y^ - l)\X{y)\^ + Ci8(y, A)^. 
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renormalization group improvement and find 

GsiiM'') = ^{M^g^.u - M^M,)Q{M^) ri{vv') (52) 

/ on 

Tt[Hb{v)Y{1 - 75)%*) (t^')] Tt[H-,.,{v')Y{1 - 1,)Hb{v)]. 

The T contribution follows from Eqs. (PB|)-(PBD replacing Vcs, Nq and Ci by 
one and the charm quark mass by the r mass. 

Table ^ shows the predictions for the total rates of these two ansatze for 
the wrong charm piece. For the numerical calculation, we used |Vc6| = 0.04. 



6 Comparison with Data 

Data on one-particle inclusive B decays are available from T(4S') machines 
and also from the CERN electron positron collider LEP. At the T(4S'), the 

total production rates of Z)(*)+ and D^*^~ or D^*^^ and Z}(*)° are measured on 
the resonance, from which one can only deduce the rates for B admixture 



(53) 



Table |I] lists the available data. In the fourth row of the table we list the 
CLEO measurement of wrong charm decays B — > DX which is again the 
average over B~^ and B^. 

Although we expect our method to work best near the nonrecoil point 
vv' = 1, we integrate the spectra in order to obtain total rates, assuming 
we can still get some insight into the bulk features. Since we calculate the 
right and wrong charm contributions separately, we have to add them in 
order to obtain the number of D mesons produced in B meson decays. This 
procedure corresponds to the usual multiplicity definition of the branching 
ratio fn], p. 570]. 

The first observable one can study is the total number of D mesons 
in B decays. This quantity is considerably smaller than the usually per- 
formed charm counting since it does not include the production of Dg mesons, 
charmed baryons and charmonium, the sum of which amounts to about 20% 
of the total charm production in B decays. Using the saturation assumption 
in Eq. (|47|), our model yields a total D rate about 20% below the experimen- 



tal result. However, Eq. (|47|) was taken from the analysis of the semileptonic 
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Mode 



Model 1 



Model 2 



Experiment 



B 
B 
B 
B 
B 



DX 

'L)*X 

DdirX 

DX 

D*X 



68.8% 
51.8% 
17.1% 
7.9% 
5.0% 



68.8% 
52.7% 
16.1% 
7.9% 
5.9% 



(87.2 ± 3.5)% 
(48.7 ±3.1)% 
(38.5 ±4.7)% 
(7.9 ±2.2)% 



& + D± 
'jj*o ^ j^*± 

D - D* 
M (input) 



B 

B 
B 
B 
B 
B 



&X 

D^X 

D*^X 

D*^X 

D%j.X 



52.0% 
16.8% 
25.9% 
25.9% 
8.5% 
8.5% 



52.3% 
16.5% 
26.3% 
26.3% 
8.1% 
8.1% 



(63.1 ±2.9)% 
(24.1 ± 1.9)% 
(26.0 ± 2.7)% 
(22.7 ± 1.6)% 
(21.7 ±4.1)% 
(16.8 ±2.9)% 



0,0 



301 



p, p, PI 

'D°-D*U).68D*^ 
- 0.32D*± 



B 
B 
B 
B 
B 
B 
B 
B 



D-£+ueX 

D^e+UfX 

D*-i+ueX 

D*^£+u,X 

D-T+UrX 

ifr+VrX 

D*\+UrX 



2.0% 
6.5% 
3.3% 
3.3% 
0.6% 
2.0% 
1.0% 
1.0% 



2.0% 
6.5% 
3.3% 
3.3% 
0.6% 
2.0% 
1.0% 
1.0% 



(2.7 ±0.8)% 
(7.0 ± 1.4)% 
(2.8 ±0.4)% 

(3.2 ±0.7)% 



0,0 



Q,0 
O'i 



B'- 


.Tfx 


B^- 


^ D-X 


B+ - 


^d'x 


B+ - 


D-X 


5°- 


^D^X 


5°- 


^ D+X 


fi+ - 


-^D^X 


B+ - 


D+X 


B^- 


.D*'X 


B^- 


^ D*-X 


B+ - 


-.D*'X 


B+ - 


D*-X 




^ D*^X 




-> D*+X 


fi+ - 


D*°X 


fi+ - 


D*+X 



31.8% 
29.1% 
60.9% 
0% 
5.7% 
2.2% 
5.7% 
2.2% 

0% 
46.8% 
46.8% 
0% 
2.5% 
2.5% 
2.5% 
2.5% 



31.8% 
29.1% 
60.9% 
0% 
6.0% 
1.9% 
6.0% 
1.9% 
0% 
46.8% 
46.8% 
0% 
3.0% 
3.0% 
3.0% 
3.0% 



Table 1: Comparison of our results with data. Model 1 uses Eq. 
uses Eq. (B3). Branching ratios are computed using tb+ = r^o 



model 2 
1.55 ps. 
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Bt(B D*X) 

Tg = — 


All 

channels 


Semi- 
leptonic 


Naive spin counting 
Model 1 [see Eq. (37)] 
Model 2 [see Eq. (|38|)] 


3 

3.04 
3.26 


3 

3.36 


Fit all data using D^*'^^, 
Data using D* and -Ddir only 
Data using Z)° and only 


1.39±0.27 
1.26±0.22 
1.92±0.46 


1.58±0.70 
1.62±0.91 
1.87±1.74 



Table 2: Predicted and measured spin counting ratio r^. For the last line 
and for the fit, the relation ( |^ ) between charge and spin counting has been 
used as a constraint. For the column "all channels," the error of the fit has 
been scaled by \Jx^ /Nf, where = 4.2 and Nf = 2. 

decays where the D and D* exclusive final states saturate only (71 ± 13)% 
of the inclusive rate B —>■ Dl^vpX. Thus Eq. (|47| ) should be replaced by 

^{vv') = \^{vv')\'^ + 67]{vv'), (54) 

where Sr] accounts for the remaining contributions. Because we use factor- 
ization, this would enhance the nonleptonic channels by a similar amount 
and hence improve the prediction for the total number of D mesons. 

Keeping this in mind, we still stick to the saturation assumption in 
Eq. ( |47| ) since 61] is unknown and the naive use of Eq. ( ^7D reproduces the 
data on B ^ D*X decays. This could be accidental, as the real problem is 
spin counting. The ratio of vector to direct pseudoscalar mesons defined in 
Eq. ( PD is expected to be about three, but experiment yields a ratio barely 
above one (see Table 0). Even in the semileptonic case, the experimental 
value is slightly lower than expected, although errors are large. Since the 
rates of pseudoscalar and vector D mesons are connected by heavy quark 
spin symmetry, this large discrepancy is difficult to understand. 

Although the effect is less pronounced than the spin counting problem, 
neutral D mesons tend to occur slightly more often than expected, see Ta- 
ble Assuming the measured value for the spin counting ratio r^, the 
measured charge counting ratios tq are generally one or two standard devia- 
tions above expectations. If confirmed, this kind of effect would point to an 
additional contribution which could for instance arise from a nonfactorizing 
topology as depicted in Fig. ^ Nevertheless, current data are still consis- 
tent with the relation (|46|) between charge and spin counting, see the fit in 
Table |. 
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Channel 


Model 


Naive 


1 


2 


Data 


Bt{B - 


> '&*°X)/Bi{B D*^X) 


1 


1 


1 


1.15±0.14 


Bt{B 
Br{B 


&tX)/Br(5 ^ Di^X) 
D^X) /Bi{B D^X) 


1 
3 


1 

3.09 


1 

3.17 


1.29±0.17 
2.62±0.24 


Bt{B - 


y &HuX) /Br{B D*^£uX) 


1 


1 


1.14±0.30 


Bt(B ^ 


Dl/yX)/Bv{B D^/uX) 


1 


1 


1.05±0.89 


Bv{B 


'dU'jyX)/Bv{B D^iuX) 


3 


3. 


20 


2.59±0.93 



Table 3: Predicted and measured charge counting ratios tq. Using the mea- 
sured value rs ~ 3/2, naive charge counting yields tq — 7/3 instead of 3. 




Figure 3: Nonfactorizing isospin- violating topology, supplying a possible ex- 
planation for an enhanced D rate. 
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Another piece of information for these decays are the D momentum spec- 



tra. These have been measured by ARGUS |jT3[ and CLEO |jTj|, we use the 
data from CLEO which is more recent and more precise. The spectra are 
momentum distributions in the rest frame of the T(4S'). However, the effect 
of the motion of the B mesons produces only a neghgible smearing of these 
spectra, and we can safely ignore this effect here. In Fig. ^ we compare the 
data obtained by CLEO |]14| with our theoretical prediction. 

The spectra show that there is indeed a problem with the ratio of pseu- 
doscalar to vector mesons. Using the saturation assumption in Eq. (^Tf ), the 
spectra of the vector mesons are described within experimental uncertain- 
ties and the problem appears with the low momentum region of the decay 
spectra for the pseudoscalar mesons. Our theoretical ansatz, especially the 
SDE, should work best in this region of small D momentum. In addition, 
the shape of the spectra in this region is mainly determined by phase space, 
which yields a behavior proportional to for small pn and constant matrix 
element. The steep rise of the momentum spectra for the pseudoscalar D 
mesons is thus difficult to understand. An investigation of slow direct pseu- 
doscalar D mesons, i.e., those that do not originate from intermediate D* 
vector mesons, would be desirable. 

This problem appears to be limited to nonleptonic decays. In Fig. ^ 
we compare the semileptonic D meson momentum spectrum measured by 
CLEO ||lOl with the theoretical predictions from [|| as well as with the sum 
of the first six contributing exclusive channels |[T^ and find no significant 
disagreement. 

There is also a measurement of the wrong charm D spectrum by CLEO 
|p!0| . Unfortunately only four bins could be measured having still substantial 
uncertainties. In Fig. ^ we see that the parton-model inspired model 2 fits 
the data better than model 1, although evidence is not yet conclusive due to 
the quality of the data. 



7 Conclusions 

In this paper we have developed a QCD based description of one-particle 
inclusive decays of the type B —>■ DX and B DX. The method we suggest 
is based on the large- A^^^ limit of QCD, allowing us to factorize certain matrix 
elements. Once factorization has been performed, one can identify pieces in 
the rates which can be treated by a short distance expansion, assuming the 
bottom and the charm quarks to be heavy. This yields a series in inverse 
powers of the parameter mfj — mc- The numerator of the expansion parameter 
is a typical QCD scale for the light degrees of freedom. Thus the expansion 



16 




Figure 4: Momentum spectra l/T^otdV / dp£, of -B — > DX in comparison with 
theoretical predictions. The sohd hne is model 2, the dashed one model 1. 
The columns refer io D^ D* and direct D mesons, the rows to a charge sum, 
neutral and charged mesons. 
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Figure 5: One-particle inclusive semileptonic D momentum spectrum in i? ^ 
Dl^viX measured by CLEO [10|. The solid line is the prediction of |Q; the 
dashed lines are sums of predictions for exclusive channels following |I7|. 



18 



o 



0.08 - 



0.05 - 



0.02 




-J I I I I I I I I I I I I I I U_J L 



2 PD 2.5 
GeV 



Figure 6: One-particle inclusive wrong charm D momentum spectrum in 



B DX measured by CLEO |jTO[ in comparison with theoretical predictions. 
The solid line is model 2; the dashed one model 1. 
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parameters are Aqcn/i^b — ^c), and as{mc) and hence corrections to 

our calculation could be fairly large, in the worst case of the order of 30%. 

We have studied the leading term of this expansion which still contains a 
number of unknown nonperturbative functions. These functions have to be 
parametrized. In the same way as in the semileptonic case, we reduce the 
number of functions appearing in the right charm contributions to a single 
nonperturbative function which can be related to the Isgur Wise function, 
once we assume that most of the rate is saturated by the decays into the two 
ground state mesons S and D*. For the wrong charm case we suggest two 
models corresponding to two different ways of contracting the spinor indices. 
Both models have a single nonperturbative form factor, which we adjust to 
the experimental wrong charm yield. 

Although our method should work best close to the nonrecoil point, we 
also calculate total rates in order to discuss the total number of D mesons in 
B decays and the spin and charge counting. While the well-known problem 
of spin counting, see, e.g., the review [|r^, is not solved by our ansatz, charge 
counting seems to work well, once we properly take the D* D decays into 
account. 

The shapes of the decay spectra into vector mesons are already described 
quite satisfactorily, in particular in the low momentum region, while the D 
meson spectra are off in this region. Furthermore, our model reproduces the 
normalization of the D* spectra, such that the spin counting problem mani- 
fests itself in a deficit of D mesons. Since the experimental rates are above 
the theoretical ones, it would be interesting to investigate which exclusive 
channels contribute in the small momentum region. 
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